u=[0,1i];
subplot(3,3,1);
plot(u);
axis([-0.5,0.5,0,1]);
for k=2:8
    m=u/3;
    uu=[m,1i/3+m*(sqrt(3)/2+1i/2),m+1i/3,2i/3+m*(sqrt(3)/2-1i/2),m+2i/3];
    subplot(3,3,k);
    plot(uu);
    axis([-0.5,0.5,0,1]);
    u=uu;
end